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Abstract 

In this paper we give, for the first time, a complete description of the dynamics of 
tilted spatially homogeneous cosmologies of Bianchi type II. The source is assumed to 
be a perfect fluid with equation of state p = (7 — where 7 is a constant. We show 
that unless the perfect fluid is stiff, the tilt destabilizes the Kasner solutions, leading to 
a Mixmaster-like initial singularity, with the tilt being dynamically significant. At late 
times the tilt becomes dynamically negligible unless the equation of state parameter 
satisfies 7 > ^ . We also find that the tilt does not destabilize the flat FL model, with 
the result that the presence of tilt increases the likelihood of intermediate isotropization. 

Key words: tilted spatially homogeneous cosmologies 

1 Introduction 

Spatially homogeneous (SH) cosmologies, that is, cosmological solutions of the Einstein field 

equations that admit a local group of isometries acting on spacclikc hypcrsurfaccs, are of 
considerable importance in theoretical cosmology and have been much studied since the 
1960s. These models can be used to analyze aspects of the physical Universe which pertain 
to or which may be affected by anisotropy in the rate of expansion, for example, the cosmic 
microwave background radiation, nucleosynthesis in the early Universe, and the question of 
the isotropization of the universe itself (see, for example, [1]). 

Spatially homogeneous cosmologies also play an important role in attempts to understand 
the structure and properties of the space of all cosmological solutions of the Einstein field 
equations, since they arc part of a symmetry-based hierarchy of cosmological models of 
increasing complexity, starting with the familiar Priedmann-Lemaitre models: 



i) 


Priedmann-Lemaitre cosmologies 


ii) 


non-tilted SH cosmologies 


iii) 


tilted SH cosmologies 


iv) 


G2 cosmologies 


v) 


Gi cosmologies 


vi) 


generic cosmologies 



The terminology used in this hierarchy has the following meaning. A SH cosmology is said 
to be tilted if the fluid velocity vector is not orthogonal to the group orbits, otherwise the 
model is said to be non-tilted [2]. A G2 cosmology admits a local two-parameter Abelian 
group of isometries with spacelike orbits, permitting one degree of freedom as regards spatial 
inhomogeneity, while a Gi cosmology admits one spacelike Killing vector field. 

An important mathematical link between the various classes in the hierarchy is provided 
by the idea of representing the evolution using a state space. The physical state of a cos- 
mological model at an instant of time is represented by a point in the state space, which 
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is finite dimensional for classes i)-iii) and infinite dimensional otherwise. The Einstein field 
equations are formulated as first order evolution equations, and the evolution of a cosmolog- 
ical model is represented by an orbit (i.e., a solution curve) of the evolution equations in the 
state space. The state space of a particular class in the hierarchy is contained in the state 
spaces of the more general classes, which implies that the particular models are represented 
as special cases of the more general models. This structure opens the possibility that the 
evolution of a model in one class may be approximated, over some time interval, by a model 
in a more special class. 

The models in each level of the hierarchy can be classified according to generality in 
various ways. For our purposes the most important of these is the algebraic classification of 
the isometry group of the SH models, the so-called Bianchi classificatior^ (see, for example, 
[3], page 112). There is also a classification of the G2 cosmologies, determined by the action 
of the isometry group, that is relevant for this paper. We refer to [4] (Table 12.4 on page 
268) for details of how this classification relates to the Bianchi classification. 

In this paper we assume that the matter content of the universe is a perfect fluid with 
equation of state p = (7 — l)Ai, where 7 is constant, the cases 7 = 1 (dust) and 7 = | 
(radiation) being of primary interest. Considerable work has been done in analyzing SH 
models, subject to this assumption, and a detailed, although still incomplete, description 
of the non-tilted models has been obtained. We refer to [4] (Chapters 6 & 7)and [5], for 
details. Much less is known about the tilted models, but it is evident that as one moves 
through stages i)-iii) in the hierarchy there is an increase in dynamical complexity and new 
features emerge at each stage. It is plausible to assume that as one moves to levels iv)-vi), 
which contain inhomogeneous models, this trend will continue. Clearly, understanding the 
dynamics at one level of complexity is a prerequisite for understanding the dynamics at a 
higher level. It is within this framework, which has a long term goal of obtaining qualitative 
information about the evolution of spatially inhomogeneous models, that the analysis of SH 
models assumes renewed importance. 

For the class of tilted SH models, the Einstein field equations have been written as an 
autonomous DE in a number of different ways (see for example, [6], [7] and [8]). Nevertheless, 
due to the complexity of the equations, a detailed analysis of the dynamics has not been 
given except in the case of a subclass of models of Bianchi type V [9]-[ll]. Our goal in 
this paper is to give a qualitative analysis of the dynamics of the tilted SH cosmologies 
of Bianchi type II near the initial singularity and at late times, using the methods of the 
theory of dynamical systems. The Bianchi II cosmologies, while very special within the 
whole Bianchi class, nevertheless play a central role since the Bianchi II state space is part 
of the boundary of the state space for all higher Bianchi types (i.e. all types except for I and 
V). We thus expect that an analysis of the dynamics of the tilted Bianchi II class will give 
insight into the dynamics of the more general tilted Bianchi classes, while providing a lower 
bound for their dynamical complexity. 

The paper is organized as follows. In Section 2 we present the evolution and constraint 

^This classification applies to SH models which admit an isometry group G3 acting simply transitively 
on the hypersurfaces of homogeneity, and includes all SH models except for the Kantowski-Sachs models, 
which admit an isometry group G4 acting on the hypersurfaces of homogeneity, but with no simply transitive 
subgroup G3. 
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equations that arise from the Einstein field equations, and in Section 3 we describe the 
stabihty properties of the equihbrium points of the evolution equations, which provide the 
basis for the discussion of the dynamics in the singular asymptotic regime in Section 4 and 
in the late time asymptotic regime in Section 5. In Section 6 we discuss the implications of 
the results. Appendix A gives details of the derivation of the evolution equations. 

Parts of the paper, in particular Sections 2 and 3, and the Appendices, are inevitably of 
a rather mathematical nature. Sections 4 and 5 are less technical, and in these sections we 
give the physical interpretations of the results. In the Introduction and in the Discussion 
we discuss the longer term goals of the research and its potential significance in a broader 
context. 

The background material needed for this paper can be found in [4]. In Appendix A 
it is assumed that the reader is familiar with the orthonormal frame formalism of Ellis & 
MacCallum [12] (see [4], Chapter 1). In addition, familiarity with some basic concepts and 
results from the theory of dynamical systems is assumed in Sections 2 and 3 (see [4] , Chapter 
4). We use geometrized units with c = 1, 8%G = 1, and the sign conventions of [4]. 



2 Properties of the evolution equations 

In order to write the Einstein field equations in a form amenable to dynamical systems anal- 
ysis, we use the orthonormal frame formalism of Ellis & MacCallum [12]. In this formalism 
the commutation functions of the orthonormal frame are used as the gravitational field vari- 
ables, which has the advantage of leading directly to first order evolution equations for the 
gravitational field. The first step is to choose the orthonormal frame to be invariant under 
the group of isometries, which implies that the commutation functions depend only on a pre- 
ferred time variable t. The second step is to choose the timelike frame vector to be equal 
to the unit normal n of the group orbits, which is thus tangent to an irrotational congruence 
of geodesies. The third and key step is to make the commutation functions dimensionless by 
dividing them by the rate of expansion of the normal congruence, which leads to variables 
that remain bounded throughout the evolution of the models. 

As described in Appendix A, the class of tilted Bianchi II cosmologies can be described 
by the following set of expansion-normalized variables: 

x= (S+,S_,Si,S3,iVi,^3), (2.1) 
subject to one constraint of the form 

g{^) = h{v3)n - SgA^i = 0, (2.2) 
where h{v3) is given by (|A.19|) , i.e. 

A(^3) = G = 1 + - 

Here S_|_, E_, Si, S3 are shear variables, Ni is a spatial curvature variable and ^3 is a tilt 
variable. These variables determine the density parameter Q and the shear parameter S 
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according to (|A.17| ), (|A.38| ) and (|A.39|) . As mentioned at the end of Appendix A, we regard 
the off-diagonal shear variables S3 and Si as representing the two tilt degrees of freedom. 

The evolution equations for the variables ( |2.1| ), as derived in Appendix A, are given 
below: 



s'_ 



-(2 - g)S+ - 3S2 + lAT^ + ^SgA^i^is 
-(2 - g)S_ + 2v^S2 - VSEj - 
-(2-g + 2y3S_)Si 

-(2-g-3S+- V3S_)S3 
(g-4S+)A^i 
v,{l - vl) 



(2.3) 



(7 - 



r(37-4-S+ + V3S_) 



where 



g = 2 (1 - l^Nl) - \G-^n [3(2 - 7)(1 - v] 
The auxiliary equation for Vl' is 

VL' = G-^[2Gq - (37 - 2) - (2 - -f)vl - 7(S+ - V3S_)^;|]f]. 



(2.4) 



The state space is the subset of defined by the constraint (|2.2|) and the inequality 
n>0, which by ([AItD , ( [OsD and (|A]3|), is equivalent to 



si + s^ 



s? + s^ + 3^iVi' < 1. 



(2.5) 



This restriction, and the fact that < 1, implies that the state space is bounded. 
The evolution equations (|2.3| ) are invariant under the transformations 



(S+, S_, Si, S3, A^i, V,) ^ (S+, S_, ±Si, ±S3, ±N,, ±vs), 

provided that the product t'3S3A'^i does not change sign. In addition, it follows that A^i, Si, S3 
and f3 cannot change sign along an orbit. Thus, without loss of generality, we can assume 
that 

A^i > 0, Si > 0, S3 > 0, and V3 > 0. 

Taking these restrictions into account, the state space T> of the tilted perfect fluid Bianchi 
II cosmologies is defined by the inequalities 



A^i > 0, n>0, <V3<1, S3 > 0, Si > 0. 

The boundary &D is obtained by successively replacing the strict inequalities in 
equalities. 

The evolution equations (|2.3|) are an autonomous DE in of the form 



(2.6) 

by 



ffx) 
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where the function f : on the right side is a rational function (note the function 

G in the denominator in q, and the form of the fg equation). For values of 7 that satisfy 
< 7 < 2, the two functions in the denominator are strictly positive on the physical state 
space V and on its boundary dT>. The DE ( p.3|) is thus smooth, indeed analytic, on the 
set V U dT>. Since this set is compact and invariant, the solutions of the DE ( |2.3|) can be 
extended for all r G M. 

The constraint ( p^.2D entails a consistency requirement, namely that the equation 5* = 
should define an invariant set of the evolution equations (|2.3|) . A straightforward calculation 
using shows that 

^' = (2g-2-S+ + y3S_)(7. (2.7) 

It thus follows (see [4], Proposition 4.1 on page 29) that g = does indeed define an invariant 
set. 

In analyzing a class of Bianchi cosmologies, one typically finds that the orbits in the 
boundary of the state space play a significant role in determining the dynamics, since orbits 
in the interior can shadow[| orbits in the boundary. For the tilted Bianchi II models, the 
boundary of the state space is the union of five disjoint invariant sets, as shown in Table 
2.1. The invariant sets i)-iv) describe familiar solutions but give multiple representations of 
them. For example, the orbits with Si = in the invariant set i) describe non-tilted Bianchi 
II cosmologies relative to a Fermi-propagated frame, while the orbits with Si > describe 
the same models, but relative to a rotating frame (see (|A.37| )). Likewise, the orbits with 
Si = S3 = in the invariant set ii) describe (nontilted) Bianchi I cosmologies relative to 
a Fermi-propagated frame, while the orbits with Si > and/or S3 > describe the same 
models, but relative to a rotating frame. 





Name 




Restrictions 


Dimension 


i) 


non-tilted non- vacuum Bianchi II 


V3 


= = S3, iVi > 0, ^] > 


4 


ii) 


non-tilted non-vacuum Bianchi I 




= = Ni, Q > 


4 


iii) 


Taub (vacuum Bianchi II) 


S3 


= = Q, Ni>0,V3<l 


4 


iv) 


Kasner (vacuum Bianchi I) 




Ni = o = n,v3<i 


4 


v) 


extreme tilt set (7 < 2) 




V3 = l 


4 



Table I: The invariant sets that comprise the boundary of the state space ( |2.6| ). 



The invariant set iii) contains the usual representation of the Taub vacuum solutions ([4] 
pages 137-8 and page 196) when Si = = V3, but it also contains multiple representations 
of these solutions relative to a non-Fermi-propagated frame. Likewise, the invariant set iv) 
contains the usual representation of the well-known Kasner vacuum solutions when Si = 
^3 = "^3 = 0, but it also contains multiple representations of these solutions. The orbits in 

■^i.e. approximate closely. See [4], page 104. 
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the remaining part of the boundary, the extreme tilt set v), do not correspond to spatially 
homogeneous cosmological solutions of the Einstein field equations since equation ( |A.21| ) 
breaks down when f 3 — >■ 1 (i.e. v^v^ — >■ 1). 

In addition to the above invariant sets, the DE (|2.3| ) has one invariant set that is not a part 
of the boundary, namely the invariant set defined by Si = 0. In view of the interpretation of 
Si, this invariant set corresponds to tilted Bianchi II models with one tilt degree of freedom. 
We shall see that the dynamics of these models is significantly simpler than the dynamics of 
the full class. 

Technical point 

It is necessary for the subsequent analysis to determine at which points the constraint 
surface ( p.2|) is singular (i.e. the points x G which satisfy g'(x) = and Vg{x.) = 0), since 
at these points one cannot use the implicit function theorem to eliminate, locally, one of the 
variables. The gradient of g is 

= (-2/iS+, -2/iS_, -2/iSi, -2/1S3 - A^i, -i/iA^i - S3, h'n) , (2.8) 



where the variables are listed in the order ([27 

It follows that the surface is singular at and only at points x given by ^3 = A^i = S3 = 
r2 = 0. Thus the surface is non-singular at all points of the state space (|2.6|), and is singular 
only on part of the boundary, i.e. a two-dimensional subset of the Kasner set. Indeed, if 
> 0, the constraint can be written in the form 



and since h'^v^) > for < 1^3 < 1 and < 7 < 2, it follows that ^3 is determined uniquely 
in terms of the other variables. 

It follows from ( p.8|) that the vectors 

e^ = J l^d^^^ .2.9) 

dxA h'{v3)Q dxA <9f 3 ' 

where A = 1,2,... ,5 and (xa) = (S+, S_, Si, S3, A''!), satisfy ■ Vg = and are linearly 
independent. These vectors thus span the tangent space to the constraint surface g = in 
M^. Equivalently, an orbit lies in the constraint surface g = if and only if its tangent vector 
is a linear combination of the e^. 

3 Equilibrium points 

In this Section we consider the local stability of the equilibrium points of the DE ( p.3|) and 
the constraint (p^.2|), i.e. points x = a that satisfy 

f(a)=0, ^(a) = 0. (3.1) 
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The local stability is determined by linearizing the DE ( p.3|) at x = a, which gives 

x' = Df(a)x, 

and finding the eigenvalues of the derivative matrix Df(a). The analysis is complicated by 
the constraint, which requires that we consider only eigenvectors that are tangent to the 
constraint surface, i.e. that are orthogonal to the gradient vector Vg{sL). We shall refer 
to eigenvalues and eigenvectors that satisfy this condition as physical. We note that if all 
the physical eigenvalues have negative (positive) real parts then the equilibrium point is a 
local sink (source), i.e. it attracts (repels) all orbits in a neighbourhood. In addition to 
isolated equilibrium points, we will also encounter arcs of equilibrium points, for which one 
eigenvalue is necessarily zero (see for example [4], Section 4.3.4). In this case the criterion for 
a local sink (source) is that all eigenvalues other than the zero one have negative (positive) 
real parts. 

We now list the equilibrium points, obtained by systematically solving equations ( |3.1| ). 
Each equilibrium point, apart from those with extreme tilt, corresponds to a self-similar 
solution of the Einstein field equations^, which we also give. 

Non-vacuum equilibrium points {Q > 0): 

i) Flat FL pomt, F 

E+ = S_ = El = S3 = iVi = t;3 = 0, 
n = l, S = 0, g = i(37-2), 0<7<2. 

Self-similar solution: the fiat FL solution. 

ii) Non-tilted point, PII 

S+ = |(37-2), iVi = |V(37-2)(2-7), S_ = Si = S3 = = 0, 
fi=A(6-7), S = |(37-2), g=l(37-2), | < 7 < 2. 
Self-similar solution: the Collins-Stewart solutions ([4], pages 131 and 189). 

iii) Tilted point, Film 

S+ = |(37-2), S_ = ^(10-77), S3 = ^av/(ll7-10)(77-10), Si = 0, 



.2 _ 3(2-7) 10 



a' = ^, f < 7 < 2. 

177 -18' ^ ' 

fi = 1^2(2172-247 + 4), S^ = 1-^2(7- 1)(97- 5), g = i(37-2) 
Self-similar solution: first given in [13]. 



*The situation is analogous to the non-tilted case ([4], Section 5.2.3). 
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iv) Line of tilted points, Clltut 

S+ = i> S- = -^' ^3 = §7^(46+ 1)(8- 36), El = ^6, 



with 



0<&<1, 7=f. 
= ^(166^-456 + 59), = -|-(26+l)(9-26), g 



Self-similar solutions: not given previously. 

v) Extreme tilted point, Pllextreme 

E+ = |, E_ = S3 = El = iVi = = 1 

O 20 v2 28 ^ 4 n ^ ^, ^ O 
^^=57' ^ =57' ?=3' 0<7<2. 

vi) Jacobs disc, J 

El = E3 = ATi = = 0, E^ + < 1 
O>0, E<1, g = 2, 7 = 2. 

Self-similar solutions: the Jacobs stiff fluid solutions ([14], page 1109). 

Vacuum equilibrium points (0 = 0): 

i) Kasner circle, K, 

E^ + E^ = 1, E3 = El = iVi = t;3 = 
f] = 0, E = l, g = 2, 0<7<2. 

Self-similar solutions: the Kasner vacuum solutions ([4], pages 132 and 188) 

ii) Kasner circle with extreme tilt, fCextreme 

E^ + E". = 1, E3 = El = iVi = 0, V3 = l 
= 0, E=l, q = 2, 0<7<2. 

iii) Kasner lines with tilt, /C^^ 

E+ = V3E_ + 37 - 4, E_ = ^[-37 + 4 ± V(37 - 2)(2 - 7)], E3 = Ei = 

Ni^O, < t;3 < 1 
= 0, E = l, g = 2, |<7<2. 

Self-similar solutions: the Kasner vacuum solutions referred to a non-Fermi-propagated 
frame. 
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If 7 satisfies | < 7 < 2, there are two Kasner lines /C^j^' which join /C to /Cextremc, while if 
7 = I or 7 = 2, the two lines coincide, with S_) = j aiid S_) = ^i, — ^ 

respectively. 



We need to know whether any of the equilibrium points are local sinks or sources. It 
turns out that for each value of 7 in the interval < 7 < 2 one of the equilibria, with > 
is a local sinkQ, as indicated in Table 3.1. The equilibria are related to one another by a series 
of bifurcations that occur as 7 varies. By inspection of the coordinates of these equilibrium 
points, we observe the following transitions: 



-7=- 1=— -7=— "1= — 

' 3 T-)TT ' 7 T-)TT ' 9 rTJ ' 9 



F -A PII U Plltnt ^ Clltilt ^ Pllextreme (3.2) 

b=0 6=1 

Range of 7 Local sink 



0<7<i 


flat FL point F 


3 ^ — 7 


non-tilted point PII 


10 < 0, < 14 
7 ^ ' ^ 9 


tilted point Plltiit 


7 = f 


line of tilted points Plltiit 


f <7<2 


extreme tilted point Pllextreme 



Table II: Local sinks in the tilted Bianchi II state space. 



We can describe the mechanisms for these bifurcations, without giving full details, as 
follows. The linearization of the evolution equation for A^i at F is 

iV( = l(37-2)iVi, 

showing that the spatial curvature variable A^^i destabilizes F at 7 = |. The associated 
eigenvector is in ( |2.9| ). The linearizations of the evolution equations for S3 and 1^3 at PII 
are 

S'3 = |(77-10)S3, t;^ = |(77-10)i;3, 

showing that and S3 destabilize PII at 7 = y. The associated physical eigenspace is 
actually one-dimensional and the eigenvector is e^. Finally, the linearization of at Plltiit 
is 

E; = |(97-14)Ei, 

showing that Si destabilizes Plltiit at 7 = ^. Stability is transferred from Plltiit to Plextreme 
through the line of equilibrium points Clltat, which exists only for 7 = ^. We shall refer to 
the bifurcation at 7 = | as the spatial curvature bifurcation, and, in view of the fact that S3 

■^With the exception of F, many of the eigenvalues are compHcated expressions in 7; their exphcit form 
is unimportant for our purposes. 
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and Si represent the tilt degrees freedom, we shall refer to the bifurcations at 7 = y and 
7 = y as the first and second tilt bifurcations. 

As regards local sources, it turns out that unless 7 = 2, none of the equilibrium points or 
equilibrium sets is a local source. This result follows from a careful analysis of the eigenvalues 
associated with the equilibrium points and sets. If 7 = 2, it turns out that a subset of the 
Jacobs disc shown as the shaded region in Figure I, is a local source. 




Figure I: The shaded region, a subset of the Jacobs disc, is a local source in the case 7 = 2. 

It is of interest to consider the invariant subset defined by Ei = 0, which describes the 
evolution of models with one tilt degree of freedom. For these models an analysis of the 
eigenvalues shows that there is an arc of the Kasner circle /C, defined by < i, S_ > 0, 
that is a local source. 

4 The late time asymptotic regime 

We have seen that for each value of 7 in the range < 7 < 2, excluding 7 = y, there is a 
unique equilibrium point that is a local sink of the evolution equations, while if 7 = ^ , there 
is an arc of equilibrium points that is a local sink. These local sinks are listed in Table 3.1, 
and the bifurcations that occur as 7 increases are displayed in equation ( |3.2|) . By definition 
of local sink, any orbit that enters a sufficiently small neighbourhood of the sink approaches 
the sink as r ^ +00. The monotone functions in Appendix B, and numerical simulations, 
provide strong evidence that for a given value of 7 the local sink is the future attractor of 
the evolution equation, i.e. all orbits, except possibly a set of measure zero, approach the 
local sink as r — +00. 

The main conclusion that can be drawn from this asymptotic result is that the dynamical 
significance of the tilt and the shear at late times increases as the equation of state parameter 
increases from to 2, as follows: 
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i) in the range < 7 < | the models isotropize, and since the deceleration parameter is 
asymptotically negative (i.e. lim q = 5(87 — 2) < 0), the models are inflationary, 

T— >+00 

ii) at 7 = I the spatial curvature destabilizes the flat FL equilibrium point, and for 7 > § 
the models no longer isotropize, 

iii) at 7 = y, the tilt destabilizes the Collins-Stewart solution, and for 7 > y, the models 
are asymptotically tilted at late times, and 

iv) if 7 > Y, the tilt is asymptotically extreme (f — » 1) at late times. 

The 7-dependent limits of the dimensionless shear scalar S and the tilt variable v, as 
defined by (|A.23|) , can be obtained from the list of equilibrium points in Section 3. In the 
physically important cases of dust (7 = 1) and radiation (7 = |), which satisfy | < 7 < y, 
the models are asymptotic to the Collins- Stewart solution i.e. they do not isotropize (S ^ 0), 
but the tilt becomes dynamically negligible {v ^ 0), at late times. 

5 The singular asymptotic regime 

As mentioned in Section 3, there is no equilibrium point or equilibrium set that is a local 
source, except in the special case 7 = 2, or unless one restricts consideration to models with 
only one tilt degree of freedom (the invariant set Si = 0). The implication of this fact is that 
a typical orbit is not past asymptotic to an equilibrium point. The situation is analogous to the 
case of non-tilted SH cosmologies of Bianchi types VIII and IX (see [4], Section 6.4), for which 
it has been shown that there exist infinite heteroclinic sequences based on the circle JC of 
Kasner equilibrium points, i.e. infinite sequences of equilibrium points on /C, joined by special 
heteroclinic orbits directed into the past (see [15] and [4], Section 6.4.2, for a geometrical 
description of these heteroclinic sequences). These heteroclinic sequences determine the 
dynamics in the singular regime (r — >■ —00) in the sense that a typical orbit shadows^ (i.e. 
is approximated by) a heteroclinic sequence as r —* —00. In physical terms, the dynamics 
of a typical cosmological model is approximated by a sequence of Kasner vacuum models as 
the singularity is approached into the past, the so-called Mixmaster oscillatory behaviour. 

The present situation is more complicated due to the existence of two Kasner circles, the 
standard Kasner circle /C, and the Kasner circle /Cextreme, with extreme tilt (i.e. 1^3 = 1). The 
heteroclinic sequences contain orbits that join two points on /C, orbits that join two points 
on /Cextreme, and orbits that join a point on /C to a point on /Cgxtreme, and vice versa. There 
are three families of orbits that join two points on the same Kasner circle. These orbits 
satisfy = 0, i.e. 

+ + T.\ + T.l + ^Nl = 1, 
have one of Si, S3 and A'^i non-zero, and have 

t'3 = or = 1. 

^In the language of dynamical systems, we are describing the a-limit set of a typical orbit. For recent 
progress in proving the existence of the a-limit set, we refer to [16] & [17. 
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The three famihes are given by 



i) N,>0, E_ = Ci(E+-2), 

ii) E3 > 0, E+ - v^E_ = C2, 
hi) El > 0, E+ = C3, 

where Ci,C2 and C3 are constants. The projections of these orbits in the E_|_E_-plane are 
shown in Figure II, with the arrows showing evolution into the past. We note that the orbits 

i) describe the Taub vacuum solutions of Bianchi type II (see [4], page 137), while the orbits 

ii) and iii) describe the Kasner vacuum solutions relative to a rotating frame (i.e. the angular 
velocity of the spatial frame is non-zero; see Appendix A). 




i) Nj >0 



KorK, 

extreme 




ii) Z3 > 



Qi 




iii) El >0 



Figure II: Vacuum orbits joining points on the Kasner circle /C or the extreme Kasner circle 
^extreme) given by + = 1. The arrows show evolution into the past. 



The orbits joining an equilibrium point on /C to an equilibrium point on /Cextreme are given 

by 

iv) E+, E_ = constant, E^ + E^ = 1, Ei = E3 = A^i = 0, 

with v'^ < (i.e. ^3 varying between and 1 for evolution into the past), or v'^ > 0. The 
direction of flow along these orbits is determined by the Kasner lines with tilt JC^, which 
depend on the value of 7, and is shown in Figure III. 
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a) Y = l b) 3^<Y<2 c) y = 2 

Figure Vacuum orbits joining points on the Kasner circle /C and the extreme Kasner circle 

^extreme- The arrows show evolution into the past. 

If < 7 < |, the flow into the past along these orbits is from /C to /Cextreme, while if 7 = 2, 
the flow into the past is from /Cextreme to /C. On the other hand, if | < 7 < 2, one family of 
orbits links /C to /Cextreme while the second family does the reverse. 

We can now describe the infinite heteroclinic sequences. In the case | < 7 < 2, the orbits 
that join successive Kasner points belong to the eight families i)-iii) with f 3 = or f 3 = 1 
and iv) with tin > or < 0. In the case < 7 < |, since lim = 1, only orbits in the 

three families i)-iii) with f 3 = 1 are permitted. 

It is not our intention to describe the detailed structure of the heteroclinic sequences, as 
has been done in the case of the non-tilted models of Bianchi types VIII and IX (see the 
references at the beginning of this section). We simply wish to point out that the evolution of 
tilted cosmologies of Bianchi type II in the singular asymptotic regime is governed by infinite 
heteroclinic sequences based on the Kasner vacuum models. 

The case 7 = 2 is exceptional in that there is a local source, namely a subset of the Jacobs 
disc (see Figure I). We thus conjecture that a typical orbit with 7 = 2 is past asymptotic to 
an equilibrium point in this set. The models described by the invariant set Si = are also 
exceptional, since there is a subset of the Kasner circle JC that is a local source (see the end 
of Section 3). We thus conjecture that a typical orbit in this set is past asymptotic to an 
equilibrium point in /C. 

We conclude this section by comparing our description of the oscillatory behaviour in the 
singular regime using dynamical systems methods in the expansion-normalized state-space 

^These invariant sets of orbits joining K, and /Cextreme first appeared in the analysis of the tihed SH models 
of Bianchi type V (see [10] and [11]). 
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with the descriptions provided by the Hamiltonian approach (see [18] and [19], pages 63-4) 
and by the so-called BKL approach ([20], pages 533-8). 

For the class of non-tilted SH models of Bianchi types VIII and IX this oscillatory be- 
haviour is described as a succession of 

i) bounces off Bianchi type II potential walls, in the Hamiltonian approach, 

ii) changes of the Kasner exponents as described by the BKL map (see [4], page 236), in 
the BKL approach, and 

iii) vacuum Bianchi II orbits (Taub orbits) linking Kasner equilibrium points in the expansion- 
normalized state-space. 

The presence of tilt leads to a new dynamical phenomenon which is due to the occurrence 
of off-diagonal shear degrees of freedom, and which is described in the above three approaches 
as follows: 

i) bounces off "centrifugal" potential walls in the Hamiltonian approach (see [21], pages 
132-4), 

ii) rotation of the Kasner axes, thereby permuting the Kasner exponents, in the BKL 
approach (see [22], pages 640-7), and 

iii) non-singular Kasner orbits that link physically equivalent Kasner equilibrium points 
(see Figure II, ii) and iii)) in the expansion- normalized state-space. 

In the expansion-normalized state space, the presence of tilt leads to an additional new 
phenomenon, namely transitions between states with zero tilt and states with extreme tilt, 
described by the non-singular Kasner orbits that link equilibrium points on K, and A^extreme 
(see Figure III). These orbits, that describe changes in the dynamical significance of the tilt 
during the oscillatory regime, do not appear to have an analogue in the other two approaches. 

6 Discussion 

In this paper we have introduced expansion-normalized variables to describe the evolution 
of tilted SH models of Bianchi type II, and have shown that with this choice of variable, the 
Einstein field equations reduce to an autonomous DE on a five dimensional compact subset 
of M^. Since the state space is compact, the orbits of the DE have a well-defined asymptotic 
behaviour as the dimensionless time r tends to ±oo. This property enabled us to give a 
detailed description of the dynamics of the above class of models in the late time regime (see 
Section 4) and on approach to the initial singularity (see Section 5). 

The description is incomplete in one respect, namely, that we have not discussed the 
kinematical properties of the cosmological fluid. For a non-tilted SH cosmology the only 
non-zero kinematical quantities of the fluid congruence are the rate of expansion and the 
rate of shear, since the acceleration and vorticity are necessarily zero [12]. On the other 
hand, in a tilted SH cosmology, all four kinematical quantities of the fluid congruence are 



15 



non-zero in general [2]. It is essential to keep in mind, however, that in our analysis, the 
dimensionless shear variables (|2.1| ) and the Hubble scalar H describe the kinematics of the 
timelike congruence that is normal to the group orbits. The kinematical quantities of the 
fluid congruence can be expressed in terms of the variables (|2.1|) by adapting the formulas 
in [2] (see sections 1 and 2). The desired relations, valid for all Bianchi types, are given in 
Appendix C, and can easily be specialized to models of Bianchi type II. 
We can draw the following conclusions from these relations: 

i) The cosmological fluid is expanding for all time. This fact follows from ( |C1| ) and 
( |C2| ), and the fact that the Hubble scalar of the normal congruence is assumed to be 
positive. 

ii) The vorticity of the fluid is zerof], but the acceleration is non-zero, provided that the 
pressure is non-zero (7 7^ 1). 

iii) If the dimensionless shear of the normal congruence is small, then the dimensionless 
shear and acceleration of the fluid congruence are small. 

It is known that the evolution of non-tilted SH cosmologies of Bianchi type II has a very 
simple description in the expansion-normalized state space: the orbit of a model joins two 
equilibrium points, or is itself an equilibrium point .0 Since equilibrium points correspond 
to self-similar solutions of the Einstein field equations, one says that the models are asymp- 
totically self-similar at the initial singularity (r —>■ —00) and at late times (r +00). The 
initial asymptotic state is typically a Kasner solution and the late time asymptotic state 
is the Collins-Stewart perfect fluid solution. We have seen that a Bianchi II model can in 
general accommodate two tilt degrees of freedom, described by the two off-diagonal shear 
variables Si and E3. Our analysis shows that models with only one tilt degree of freedom 
(i.e. Si = 0) are still asymptotically self-similar. The initial asymptotic state is again a 
Kasner state, implying that the tilt is not dynamically significant near the initial singular- 
ity. The late time asymptotic stateQ depends on the equation of state parameter 7, and is 
described by the sequence of bifurcations 

a subset of ( p.2|) . If the second tilt degree of freedom is activated (i.e. Si 7^ 0), we have 
seen that evolution is Mixmaster-like near the initial singularity, and hence the models are 
no longer asymptotically self-similar in the singular regime. The late time asymptotic state 
is described by the sequence of bifurcations (p.2[), showing that if 7 satisfies ^ < 7 < 2 the 
tilt becomes extreme. For models containing dust (7 = 1) or radiation (7 = |), however, 
the models are asymptotically self-similar, with the Collins-Stewart solution as the late time 
asymptote. 

^This result can be inferred directly from [2] (see Theorem 3.1). 

^See [4], Section 6.3.2. We refer to [7] (page 53) and [23], for alternative descriptions. 
^°This asymptotic behaviom- has previously been discussed by [13]. 
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The above results lead to some interesting comparisons concerning the ways in which var- 
ious anisotropics, in particular, spatial curvature, magnetic fields and tilt, affect the dynamics 
of SH cosmologies. 

It is known that the anisotropic spatial curvature in a SH cosmology, as described by the 
variables A^^, ot — 1,2,3, affects the stability of the circle of Kasner equilibrium points (see 
[4], page 132). If all of the are non-zero, corresponding to models of Bianchi types VIII 
and IX, all Kasner points become saddles, leading to the existence of infinite heteroclinic 
sequences into the past, and hence an oscillatory singular regime. It has also been found 
that magnetic fields affect the Kasner equilibrium points in an analogous way, although 
the mechanism is somewhat different in that a magnetic field with more than one degree 
of freedom generates off-diagonal shear components. We refer to [24]-[27] for details. The 
present paper shows that in SH models of Bianchi type II, the tilt degrees of freedom affect 
the Kasner circle in the same way by generating off-diagonal shear components, leading to 
an oscillatory singular regime. Since Bianchi II orbits are contained in the boundary of the 
state space of models of more general Bianchi type (all except types I and V), we expect 
that this mechanism will operate quite generally for tilted Bianchi models. 

There are also analogies in the way in which anisotropic spatial curvature, magnetic 
fields and tilt affect the dynamics of SH cosmologies at late times. If the equation of state 
parameter satisfies < 7 < |, then all orbits are asymptotic in the future to the fiat FL 
equilibrium point F, implying that the models isotropize. Since in this case the limit of the 
deceleration parameter q is negative ( lim q — ^(37 — 2) < 0), the models are inflationary. 

For each of the three types of anisotropy, bifurcations occur as 7 increases, the result of 
which is that these anisotropics influence the dynamics to an increasing degree at late times. 
First, the spatial curvature destabilizes the flat FL equilibrium point F at the value 7 = |- 
by the creation of new equilibrium points with non-zero spatial curvature. At larger values of 
7, a magnetic field in models of Bianchi types II and VIq destabilizes these new equilibrium 
points, with the result that the magnetic field becomes dynamically significant (see [24]- 
[25]). In this paper we have shown that the tilt in a Bianchi II model acts in a similar way, 
becoming dynamically significant if 7 > y. 

There is, however, an important difference between the effect of tilt in a SH model of 
Bianchi type II and the effect of spatial curvature and magnetic fields. We have mentioned 
that the spatial curvature destabilizes the fiat FL point F at 7 = |. It has also been shown 
that a magnetic field leads to a second destabilization of F at 7 = | (see [26]). On the other 
hand, the analysis in Section 3 shows that in a Bianchi II model the tilt does not destabilize 
F. In other words, the flat FL equilibrium point F is stable with respect to the off-diagonal 
shear degrees of freedom that arise in the presence of tilt. In the language of dynamical 
systems the tilt increases the dimension of the stable manifold of the equilibrium point F 
making it more likely that orbits will pass close to F. The physical significance of this 
result is that in Bianchi II models tilt increases the probability of intermediate isotropization 
with il ^ 1. Whether this result is true in general for tilted SH models requires further 
investigation. 
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Appendix A 



In this Appendix we derive the evolution equations ( p.3|) -( pl4D for the tilted SH models of 
Bianchi type II. We begin by giving the evolution equations for a general SH model, in terms 
of expansion-normalized variables defined relative to the timelike congruence that is normal 
to the group orbits. These equations follow directly from the standard orthonormal frame 
equations for SH models as given in [4] (see equation (1.90)-(1.100)), as we now describe. 

We introduce a group-invariant orthonormal frame {eo,eQ,}, where eo = n is the unit 
normal to the group orbits. The non-zero commutation functions are ([4], page 39) 

and the energy-momentum tensor 

Tab = firiarib + 2q(^anb) + p{gab + rianb) + -Kab (A.l) 
is described by the source terms 

relative to the chosen orthonormal frame. 

The expansion-normalized commutation functions and source terms are defined by 

V - /\r - """^ 4 - P - 9^ 

We also require the shear parameter E, defined by 

(A.4) 



3//2' 

equation (5.17)) and the Hubble-normalized spatial curvature variables, defined by 



•5.. = ^, K = -^,. (A.5) 



It follows from (pT^ ) and equations (1.28), (1.94) and (1.95) in [4] that 

= (A.6) 

and 

So.p = B^(, - - 2e^'\^Np)^A,, K = ^^S/ + A^A^ (A.7) 

where 
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We introduce the usual dimensionless time variable r according to 



dt 1 / . X 

where t is clock time along the normal congruence (see [4], page 113). In making the 
transition to expansion-normalized variables, the deceleration parameter q plays an essential 
role (see [4], pages 112-3). This dimensionless scalar determines the evolution of the Hubble 
scalar H according to 

H' = -il + q)H, (A.9) 

([4], equation (5.11)). Raychaudhuri's equation ([4], equation (1.90), in conjunction with 
( |A.3| ), ( |A.4D , ( [A.8| ) and ( |A.9| ), leads to the following algebraic expression for q: 

Q = 2S2 + i(n + 3P). (A.IO) 

The Einstein field equations ([4], equations (1.91)-(1.93)), Jacobi identities ([4], equations 
(1.96)-(1.98)) and contracted Bianchi identities ([4], equations (1.99) and (1.100)), in con- 
junction with ( |A.2|) , ( |A.3|) , (|A.5|) , ( |A.8|) and ( |A.9|) , now yield the following sets of equations. 



Gravitational evolution equations: 

= -(2 - q)^^p + 2e^\^i:(,)^R, - S^p + n,^, (A.ll) 

= QNap + 2S(/Ar^)^ + 2e''^ ^^N^),R., (A.12) 

A'^ = qA^ - J:£Afs + erA^K, (A.13) 

Source evolution equations: 

n' = {2q - i)n - 3P - |s/n/ + |a„q" (a.m) 

= 2{q - 1)Q„ - S/Q^ - e^R^Qu + 3A^n,^ + srN/Up,. (A.15) 
Algebraic equations: 

N/A^ = 0, (A.16) 

n = - K, (A.17) 

Qa = 3S„M^ - er^fNfs.. (A.18) 
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It should be noted that the source evolution equations (i.e. the contracted Bianchi identi- 
ties) are a consequence of the gravitational evolution equations and the algebraic equations, 
and hence contain no additional information. It is, however, convenient to use them as 
auxiliary equations. 

It should also be noted that equations (^]Tl])-(^]T3D, with (pTTOl ), { \KJ7\ j and { \KJEi ) 
do not form a fully determined system of evolution equations. First, there is no evolution 
equation for the variables Ra that represent the angular velocity of the spatial frame {e^}. 
One can in fact use the freedom in the choice of spatial frame, i.e. an arbitrary time-dependent 
rotation, to introduce a non-rotating spatial frame = 0). This choice is not usually the 
most convenient one, however. In addition, there is neither an evolution equation nor an 
algebraic equation for the isotropic pressure P and the anisotropic stress matrix Ha/s- Thus 
in order to obtain a fully determined system one has to specify the source and fix the spatial 
frame. 

We now consider the case of a tilted perfect fluid, with stress energy tensor 



Tab = fJ-UaUb + p{gab + UaUb) , 

and equation of state 

p= (7- 

The 4-velocity u can be written in the form 



u 



(A.19) 



(A.20) 



(A.21) 



where the spacelike vector v is orthogonal to the normal vector n, and satisfies < fftf ^ < 1. 
The vector v is called the tilt vector of the fluid. It has components (0, v"") relative to the 
orthonormal frame {n, Bq,}. 

We can now express the source terms P, Qa and Ilap in terms of the density parameter 
VL and the tilt vector Wq,, as follows. First, equations ( [A.l[ ) and ( |A.19| )-( [A.21| ) imply that 



G 



fx 



1 — v' 



where 



G = 1 + (7 - 1)^;^ v^ = Vo,v"<l. 

By using ( [A.22|) , in conjunction with (|A.3| ), we now obtain 

P=IG-' [3(j-l)(l-v') + jv'] n, 
Qa = S'yG^'^Qva, 

The algebraic equation ( A.18| ), with ( [A. 25 ), assumes the form 



(A.22) 



(A.23) 



(A.24) 
(A.25) 
(A.26) 



(A.27) 



20 



Equations ( A.17|) and ( A. 27 ) express the source terms Q and Va algebraically in terms of the 
gravitational field variable S^^, A^^^ and A^. 

It is, however, convenient to use the source evolution equations ( |A.14| ) and ( [A.15| ) to 
obtain evolution equations for f2 and Va- On substituting (|A.24|) -( |A.26D in ( |A.14| ) and ( |A.15| ) 
and rearranging we obtain 



fi' = G-^ \lGq - (37 - 2) - (2 - 7)t;2 - 7S„^i;°^;'^ + 27v4„w"] fi. 



(A.28) 



au' 



dQ 



[(37-4)(l-^;2) + (2-7)S^y^;^ 



[1 - (7 - \)v^\ 
+ {(2 - 7) - (7 - 1)(1 - v')}A,vf'\ - 



(A.29) 



Finally, we substitute (^]2|) in ([AlOD to express q in terms of the gravitational field 
variables: 



g = 2E2 + \G-^^ [(37 - 2)(1 - v'') + 27^;2] . 
We note that ([A.17| ) can be used to write q in the form 

q = 2(\-K)- iG'-il][3(2 - 7)(1 - v") + 271;' 



(A.30) 



(A.31) 



For the class of Bianchi II cosmologies we have A^ = 0, and in addition we can choose 
the spatial frame {e„} to be an eigenframe of the matrix A^^^^, with 

Nu ^ 0, Nap = otherwise. (A.32) 

These restrictions, in conjunction with the constraint ([A.27|) , imply 

vi = 0, 

assuming f2 > and 7 > 0. In addition, equation ( |A.12[ ) gives R2 = — S13, and R3 = S12. 
We are free to perform a rotation in the 23-plane to get[^ 



V2 = 0, V3^ 0. 



The constraint ( |A.27| ) now yields 



Si3 = = i?2, S12 ^ 0, 



(A.33) 



Using these results the E'^g equation implies 



Ri 



^23, 



^^For details of the derivation of this equation we refer to [28], equations (65) and (94). 
^^This choice is the one made in [2], page 223. 
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and at this stage Ra is uniquely determined in terms ofEap- 
We now relabel the variables as follows: 

^+ = 1(^22 + S33), S_ = ^(S22 — S33), 

^1 = 73^23, S3 = ^Si2, 

iVi = iVn. 

The set of independent expansion-normalized variables is 

(S+,S_,Si,S3, A^i,i)3), 
subject to one constraint ( [A.33| ), which we now write in the form 

h{vs)n = S3iVi, (A.34) 



where 



and G is given by 



M"3) = (A.35) 



a = ! + {-,- l)v^, (A.36) 



as follows from ([A.23|) and the restrictions on Va- 

It should be noted that the off-diagonal shear components Si and S3 determine the 
angular velocity of the spatial frame according to 

(i?„) = v^(Si,0,S3). (A.37) 

The shear parameter has the simple form 

= + S2 + + ^l, (A.38) 

and the curvature parameter K is given by 

K = ^Nl (A.39) 

as follows from ( |A.7| ) and ( |A.32| ). The deceleration parameter q, as given by ( |A.31| ), now 
assumes the form 

g = 2 (1 - ^iV,^) - iG-'n [3(2 - 7)(1 - ^^3') + 27^1] • (A.40) 

The evolution equations (|2.3|) -( p^ are now obtained by specializing equations ( |A.11| ), ( |A.12| ), 
( |A.28| ) and (|A.29|) using the restrictions obtained above. 

We conclude this Appendix with an important remark concerning the interpretation of 
the expansion-normalized variables. In Bianchi II models there are only two tilt degrees of 
freedom instead of the customary three (see [4], Table 9.5, page 211). The reason is that 
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the matrix N^p has two zero eigenvalues which, in conjunction with the constraint ( A. 27 ), 
imphes that 

N^^vp = 0, 

i.e. the tilt vector lies in the two-dimensional null eigenspace of A^^^^^. Our choice of frame, 
which leads to fi = f2 = 0, obscures the fact that there are two tilt degrees of freedom. 
The constraint (|A.34|) shows that the shear variable S3 represents one of the tilt degrees of 



freedom. It transpires that the shear variable Si represents the second tilt degree of freedom, 
since if Si = 0, the tilt vector is an eigenvector of S^^, i.e. 



'la^p^Vf, = 0, 



which fixes its direction uniquely. 



Appendix B 

In this Appendix we give some functions that are monotone along the orbits of the evolution 
equations, depending on the value of 7. 



i) Consider the function 



fi3 ' 
where 

a- 



(1_ ^2)1(2-7)' 

and G is given by equation (|A.36|) . It follow^ from the evolution equations (|2.3| )-(p^ 

that 

Z' = 3(97 - U)Z. (B.2) 

It thus follows that if 7 satisfies 0<7<2, 77^^^ then Z is a monotone function 
in the invariant set defined by > 0, > 0, Si > and S3 > 0. Furthermore, if 
< 7 < ^, then 



lim V3S1S3 = 0, lim (1 - vl)n = 0, (B.3) 

T — ^+00 T — ^ — 00 



and if ^ < 7 < 2, then 



lim (1 - v^)n = 0, lim ^;3SiS3 = 0. (B.4) 



^■^The calculation is simpler if one uses the constraint (2.2) to express f2 in terms of the other variables. 
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ii) For some values of 7 the tilt variable is monotone. If 7 satisfies < 7 < |, the 

factor 37 — 4 — S+ + \/3S_ in the DE for t>3 is negative (see ( p.3| )), and hence is 

2 

3' 

lim V3 = 0, lim ^3 = 1. (B.5) 



monotone decreasing into the future. It thus follows that if < 7 < ^ 



If 7 = 2, the DE for ^^3 simplifies to 

v's = {2-E+ + v^S_)t;3. (B.6) 

Since 2 — S+ + v^S_ > unless S_) = — ^ j , it follows that V3 is monotone 
increasing into the future. Thus, if 7 = 2, 

lim V3 = 0. (B.7) 

T—^ — 00 

The DE (|B.6|) also implies that if 7 = 2, the extreme tilt set given by t'3 = 1 is not an 
invariant set, with the result that orbits can pass from the region ^3 < 1 to the region 
f 3 > 1 in state space. 

iii) It follows from (g^Sj) and (|2]|) that 

(/?n)'=[2g-(37-2)](/?n), (B.8) 

where 

The expression (|A.31|) for q can be rearranged to give 

2q - (37 - 2) = 3(2 - 7)S^ + (2 - 37)^ + 7(4 - 3^)G-^v^ - n, 

where K is given by ( |A.39| ). Thus, if 7 satisfies < 7 < |, PQ is monotone increasing, 
while if 7 = 2, PQ is monotone decreasing. It follows that if < 7 < |, then 

lim n = l, lim = 0. (B.9) 



~oo 



We note that the limits (§J), (|R4D , (§J), and (|R9D provide support for the 

claims made in Sections 3 and 4 concerning the asymptotic behaviour. 



Appendix C 

In this appendix, we give the equations relating the kinematical quantities of the fluid con- 
gruence in terms of the dimensionless commutation functions ( |A.2| ) and ( |A.3| ) associated 



with the normal congruence. The fluid kinematical quantities, namely, the acceleration 
the vorticity vector Ua and the rate of shear tensor cr^h^'^, are normalized using the Hubble 
scalar -fffluid of the fluid congruence: 



Un COn 



^fluid 



— T7 ; — 77 , '^ab " 77 • 

-n fluid -"fluid -"fluid 
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We now give the desired relations, which are obtained by adapting the results of [2] (see 

between our til 

Va = (tanh/?)ca, 
tanh/3 = V. 



Sections 1 and 2)R. Note the relation between our tilt vector Va and their vector Ca- 



where 



Hubble scalar. 



where 



-f^fluid = BH, {CI] 



Acceleration: 



Ua = — 1) cosh /3 fa, Uq = —v"Ua 



Vorticity: 



= l^iNa % + e^^-'v.A, + cosh^ pm'^v.v^v^), Wo = -v'^W^ 

ZlD 



Shear. 



^aP + = — + ^afs) - COsh^ (4 - 3'y)VaVf3 

cosh (3 



yfi _ _yfi yf^ _ yf^ ,,«,,/3 

^Oa — ^aP^ ' ^00 — ^aP^ ^ ■ 
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